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A new theory of coherent structure in wall turbulence is presented. The theory is the 
first to predict packets of hairpin vortices and other structure in turbulence, and their 
dynamics, based on an analysis of the Navier-Stokes equations, under an assumption of 
a turbulent mean profile. The assumption of the turbulent mean acts as a restriction on 
the class of possible structures. It is shown that the coherent structure is a manifestation 
of essentially low-dimensional flow dynamics, arising from a critical layer mechanism. 
Using the decomposition presented in McKeon & Sharma (J. Fluid Mech, 658, 2010), 
complex coherent structure is recreated from minimal superpositions of response modes 
predicted by the analysis, which take the form of radially-varying travelling waves. The 
leading modes effectively constitute a low-dimensional description of the turbulent flow 
which is optimal in the sense of describing the resonant effects around the critical layer 
and which minimally predicts all types of structure. The approach is general for the full 
range of scales. By way of example, simple combinations of these modes are offered that 
predicts hairpins and modulated hairpin packets. The example combinations are cho- 
sen to represent observed structure, consistent with the nonlinear triadic interaction for 
wavenumbers that is required for self-interaction of structures. The combination of the 
three leading response modes at streamwise wavenumbers 6, 1, 7 and spanwise wavenum- 
bers ±6, ±6, ±12 respectively, with phase velocity 2/3, is understood to represent a 
turbulence "kernel" , which, it is proposed, constitutes a self-exciting process analogous 
to the near-wall cycle. Together, these interactions explain how the mode combinations 
may self-organise and self-sustain to produce experimentally observed structure. The 
phase interaction also predicts important skewness and correlation results known in the 
literature. It is also shown that the very large scale motions act to organise hairpin-like 
structures such that they co-locate with areas of low streamwise momentum, by a mech- 
anism of locally varying the shear profile. These energetic streamwise structures arise 
naturally from the resolvent analysis, rather than by a summation of hairpin packets. In 
addition, these packets are modulated through a "beat" effect. The relationship between 
Taylor's hypothesis and coherence is discussed and both are shown to be the consequence 
of the localisation of the response modes around the critical layer. A pleasing link is made 
to the classical laminar inviscid theory, whereby the essential mechanism underlying the 
hairpin vortex is captured by two obliquely interacting Kelvin-Stuart (cat's eye) vortices. 
Evidence for the theory is presented based on comparison to observations of structure 
reported in the experimental, transitional flow and turbulent flow numerical simulation 
literature. 
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1. Introduction 

The existence of coherent vortical structure in wall turbulence has been known for 
many decades, but the exact definition of coherent structure remains controversial, be- 
yond the requirement that coherent motion has significant correlation in space and time 
(Robinson 1991). Attempts to describe additional features of the flow solely in terms 
of the velocity fields associated with the structures that are observed have had limited 
success. One notable exception is the attached eddy model proposed by Townsend (1976) 
and subsequently developed by Perry and coworkers, e.g. (Perry & Chong 1982; Perry 
et al. 1986; Perry & Marusic 1995). One weakness of the attached eddy model is that 
it does not incorporate the dynamics of the structiires. Over time, a split has developed 
in our understanding of how wall turbulence scales with increasing Reynolds number. 
One view derives from statistical measures of the velocity field and scaling arguments 
based on dimensional analysis in physical and spectral space. A second view focuses on 
the dynamics and spatial organisation of coherent structure. In the following, we refer 
to these views as "statistical" and "structural" respectively. The analysis of McKeon 
& Sharma (2010) successfully modelled some selected structural features and Reynolds 
number scaling of velocity fluctuations in turbulent pipe flow. In this contribution, we 
use that framework to reconcile the distinct statistical and structural views, by demon- 
strating that the velocity response modes of the previous analysis are associated with 
strong evidence of coherent vortical structure. The superposition of just a few such modes 
captures the essence of even complex assemblies of individual structures such as hairpin 
vortices, and the relationship between large scale coherence and individual vortices. In 
particular, we propose that hairpin vortices, hairpin packets and amplitude modulation 
arise naturally from the linear superposition of a small number of response modes. The 
characteristics of such superpositions quickly become complex as more response modes 
are added. However, the essential processes captured by simple packets are an innate 
feature of wall turbulence and we can extract them from the governing equations. This 
permits very economical descriptions of classes of coherent structure with representative 
mode superpositions. 

The structure of the paper is as follows. Firstly, we continue this introduction with a 

brief review of the current understanding of coherent structures in wall turbiilence rele- 
vant to the present study. We will then outline related modeling efforts and the connection 
with the dynamical systems viewpoint. In section 2 we present the theoretical develop- 
ment and analysis that provide the modes on which the work is based. Section 3 presents 
mode combinations representative of important characteristics of wall turbulence. Sec- 
tion 4 is a presentation and analysis of structure arising from the mode combinations, 
in particular, hairpins and hairpin packets. Their structural organisation and amplitude 
modulation across scales is also discussed. Next, section 5 discusses the results in context 
of the current understanding of structure in wall turbulence and transitional flow. In the 
final section, conclusions are presented. 

1.1. Observations of coherent structure in wall turbulence 

The review of Smits et al. (2011) identifies and describes four different classes of coherent 
structure, or organized motion, under current investigation in the literature, namely: 

(a) the streaks and quasi-streamwise vortices associated with the near- wall cycle; 

(6) hairpin or horseshoe vortices spanning from the wall into the overlap region in the 
mean velocity and perhaps beyond; 

(c) Large Scale Motions (LSMs), often associated with packets of hairpin vortices and 
stretching several outer lengthscales in the streamwise direction; 

(cf) and Very Large Scale Motions (VLSMs) or superstructures, extending for approx- 
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imately ten outer lengthscalcs in the strcamwise direction and bearing an increasing 
proportion of streamwise energy and Reynolds shear stress for increasing Reynolds num- 
ber. 

The main focus here will be on the origin and interactions of the second class of hairpin 
vortices, understood to have the distinctive, symmetrical hairpin shape only in the sta- 
tistically representative or averaged sense, and subsequently make inferences about the 
likely relationships between assemblies of these structures and the last two long classes of 
structure. Readers with an interest in predictions regarding near-wall motions and other 
structure are referred to McKeon et al. (2010) and McKeon et al. (2013). 

The seminal, classical works of Theodorsen (1952) and Head & Bandyopadhyay (1981) 
utilized flow visualization to infer correctly many now-accepted features of coherent vor- 
tices. These include the geometrical similarity of the averaged structures, downstream- 
inclined horseshoes that reach up into the flow from the near-wall region and the dom- 
inance of hairpin heads with rotation in the sense of the mean shear and clustering of 
vortices into larger packets that could approach the outer lengthscale in cross-stream di- 
mension. These conclusions were subsequently developed through the analysis of Robin- 
son (1991) which resembled hairpins in the statistically averaged sense. Falco (1977) 
presented additional observations from smoke visualization and hot wire measurements, 
although these were interpreted somewhat differently by the author. 

A concerted effort by the group of Adrian, as summarized in Adrian (2007), has con- 
firmed the existence of hairpin vortices, and demonstrated the existence of hairpin packets 
throughout the turbulent layer, culminating in LSMs. Kim & Adrian (1999) also iden- 
tified VLSMs and hypothesized that they are made up of aligned LSMs, in a picture 
that begins with very small hairpins and ends with structures much longer than the 
thickness of the turbulent layer. The statistical characteristics of hairpin-like vortices in 
a turbulent boundary layer have been explored by Carlier & Stanislas (2005) (and sub- 
sequent work) and Gao et al. (2011), and linear stochastic estimation of likely vortical 
geometries and spatial relationships have been documented by Adrian et al. (2000 a), 
Natrajan et al. (2007), Lee & Sung (2011) and others, in both experiments and simu- 
lations of wall turbulence. Ganapathisubramani et al. (2003) confirmed the dynamical 
significance of hairpin vortices, quantified by the association of more than 25% of the 
total Reynolds stress arising from Q2 events with hairpins occupying less than 4.5% of 
space (Ganapathisubramani et al. 2003). 

Despite a wealth of observational evidence for the existence of coherent hairpin vor- 
tices, there remains some controversy over their role in fully-developed turbulence, given 
the well-known occurrence of such structure in transitional flows and their apparent ab- 
sence in the very high Reynolds number boundary layer of the atmospheric surface layer 
under near thermally neutral conditions (Morris et al. 2007). Correct identiflcation of 
vortical structure in a complex, straining and rotational background velocity fleld has 
been the focus of much attention in recent years, as summarized in the works of Jeong 
& Hussain (1995) and Chakraborty et al. (2005), who point out the strengths and weak- 
nesses of different identification criteria in a range of flows most importantly associated 
with separating out coherent rotation in a field with a background vorticity (shear) likely 
of a similar magnitude to the (connected) vorticity associated with structural features. 

A proposal for the relationship between the hairpin vortices, LSMs and VLSMs was 
given by Kim & Adrian (1999) who described a possible organization in which LSMs 
consisted of hairpin vortices aligned into packet structures, which themselves aligned 
to form very large scale structures. Many authors have subsequently confirmed char- 
acteristic aspects of this description, in particular noting the common observations of 
wall-normal vortices "straddling" regions of constant low streamwise momentum, e.g. 
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Ganapathisubramani et al. (2003). The recent work of Dennis & Nickels (20116, a) de- 
veloped these results in three dimensions, by projecting three-component time-resolved 
Particle Image Velocimetry (PIV) measurements performed in the cross-stream plane 
into the streamwise direction by use of Taylor's hypothesis. They present elegant recon- 
structions of hairpin vortex packets and their collocation around extended regions of low 
streamwise momentum, constructed using averaging of the velocity field conditional on 
various swirling strength criteria. 

Guala et al. (2006) have noted that the VLSMs contain more than half of both the total 
kinetic energy of the streamwise velocity fluctuations and the Reynolds shear stress, at 
least in the moderate Reynolds number pipe flows of their studies. The subsequent studies 
of Hutchins & Marusic (20076) in the turbulent boundary layer have identified the further 
increasing energetic importance of the VLSMs as the Reynolds number increases, while 
Mathis et al. (2009) provide a comprehensive report on the influence of the very large 
scales on the smaller, near-wall turbulence activity in the form of apparent amplitude 
modulation of the envelope of small scale streamwise turbulent fluctuations by the large 
scales. 

1.2. Models for generation, distribution and, evolution of hairpin vortiees 

Existing explanations of coherent structure in both wall turbulence and transition may 
be classified as either physically-argued models, those based on linear amplification of 
some kind, or those arguing from the association with exact, self-sustaining solutions to 
the Navier-Stokes equations. 

As mentioned above, models for the origin, dynamics and evolution of coherent struc- 
ture are generally sparse in the literature, with the notable exception of the near-wall, 
self-sustaining cycle elucidated by Waleffe and co-workers, e.g. Waleffc (1997). Signifi- 
cant progress has also been made in terms of transient growth analyses, e.g. del Alamo 
& Jimenez (2006); Cossu et al. (2009). 

In terms of the hairpin vortices proposed by Theodorsen (1952), the attached eddy 
model originated by Townsend (1976), the dynamical model of Smith et al. (1991) and 
the hairpin vortex/packet paradigm developed by Adrian (2007) all address some aspects 
of their origin and development, but a complete picture appears to elude the community. 
We note also here the potential of the hierarchies of scales present in the mean momen- 
tum balance analysis of Klewicki et al. (2007), which have characteristics reminiscent of 
Townsend's attached eddies. 

An important conceptual contribution to our understanding of the origin of struc- 
ture has been made by workers looking to find self-sustaining solutions to the Navier- 
Stokes equations, either at asymptotically high Reynolds number, or "edge states" at the 
laminar-turbulent boundary in transitional fiow. It seems likely that solutions existing 
in transitional flow continue in some way to turbiilent Reynolds numbers. In this view, 
the phase-space flow around some of these solutions forms an attractor, concentrating 
the flow into regions that correspond to the coherent structures observed in turbulent 
flow. Thus the sohitions shepherd any ensemble of flow trajectories to give the distinc- 
tive statistical and structural characteristics of turbulence. The work in this area has 
concentrated on finding these exact solutions and has associated these with archetypal 
structure. To efficiently and effectively enumerate the concentrations of flow trajectories 
in the region of state space near these solutions is to "sketch" turbulence. There is some 
experimental evidence that this understanding is essentially correct in at least the case of 
transitional pipe flow (MuUin 2011). These solutions may be found by various methods, 
both numerical and analytical. It should be noted that the numerical approaches in- 
volve an iterative search and thus can be computationally expensive. It is hoped that the 
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present analysis will inform and facilitate the search for similar self-sustaining solutions 
at higher Reynolds numbers than has been possible up to now. 

In a pipe, solutions have been found by Wcdin & Kerswell (2004); Kerswell (2005); 
Duguet et al. (2008, 2010); Pringle et al. (2009); these have tended to take the form 
of self-interacting rolls and streaks. A review of such work is given by Eckhardt et al. 
(2007). Kerswell (2005) and others find that the travelling wave solutions are stationary 
states when viewed from the reference frame moving with the phase speed of the wave 
structure. That such solutions are time quasi-periodic suggests that a frequency-domain 
treatment such as that applied in McKeon & Sharma (2010) is appropriate. Localised 
solutions may have a more complicated spectral signature. 

The theoretical work of Benncy & Bergeron (1969); Smith & Bodonyi (1982) and 
others (see the review of Maslowe (1986)) found neutrally stable (self-exciting) nonlinear 
solutions for the; inviscid (equations. These asymptotic nonlinear solutions have analogiies 
at lower Reynolds numbers, for example the vortex- wave interactions found numerically 
by Hall & Sherwin (2010), and the near- wall cycle (Waleffe 1997, 2003). 

Most relevant to the current work, exact solutions appearing to consist of interacting 
hairpins and streaks were recently observed in various flows by Itano & Generalis (2009); 
Generalis & Itano (2010); Gibson et al. (2009); Deguchi & Nagata (2010) and Cherubini 
et al. (2011). 



2. Modelling and analysis: Turbulence in the frequency-domain 

In the following, we develop a formulation of the Navier-Stokes equations designed to 
examine the highly selective response of turbulent pipe flow to (internal) forcing. The flrst 
part of the analysis is a short recapitulation of the model presented by McKeon & Sharma 
(2010) with some improvements in notation. The second part makes some new observa- 
tions about the interaction of the nonlincarity with the model and discusses conseqiiences 
and some possible extensions of the development. The conceptiial picture underlying the 
analysis is one of a nonlinearity "feeding back" to excite the linear dynamics, which in 
turn drive the nonlinearity. As siich, we concentrate on the linear amplification aspect of 
the turbulent process with the understanding that the velocity fields associated with indi- 
vidual response modes may be linearly superposed. The idea is that the linear dynamics 
in the nonlinear feedback loop are so selective in their response as to bc^ tlu^ primary 
mechanism responsible for selecting structure in wall turbulence. To provide closure, we 
assume a turbulent mean velocity profile, thus restricting the search for structures to 
those that are supported by the assumed profile. By considering only velocity fields that 
are periodic in time and space, the search is also restricted to structures that do not 
decay and so, in turn, will contribute to the mean profile. Thus, the chain of reasoning 
begins with the observed mean and infers flow structures that are consistent with it. We 
find that the linear mechanism within the nonlinear feedback loop of the turbulent flow 
is highly selective (in a sense we quantify), with the consequence that analysis of a linear 
operator yields structures that are similar to those observed in wall turbulence. 

2.1. Formulation of the problem 

The non-dimensional Navier-Stokes equations for incompressible pipe flow with constant 
viscosity are given by 

ata = -Vp-u- Vu-F^V^u, V-u = (2.1) 
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Figure 1. A schematic of pipe geometry and nomenclature. 



For convenience of comparison across different studies, these equations of motion are 
non-dimensionalised witli respect to the pipe diameter and the bulk, volume- averaged 
velocity, Ubuik, so that the Reynolds number in Equation 2.1 is 



with D = 2R the pipe diameter, p the pressure field and u the kinematic viscosity. 

We retain the boundary layer terminology by fixing y = 1 — r, and u, v{= — u') and 
w as corresponding to the streamwise, wall-normal and azimuthal velocities such that 
u = {v',w,u), as shown in figure 1. 

Next, the problem is Fourier transformed in the streamwise, azimuthal and tempo- 
ral directions. The frequency-domain approach permits computational searches at high 
Reynolds number without neglecting viscosity and avoids the question of whether initial 
conditions arc "forgotten" . Wc use a different non-dimcnsionalisation when considering 
the critical layer effects because we find comparison to the centreline velocity (rather 
than bulk velocity) to be easier to work with in this context. Accordingly we non- 
dimcnsionalise variables with the pipe radius and turbulent mean ccntcrlinc velocity, 
such that k = k'R, n = n'R and ui = uj'R/Ucl, where the prime denotes the dimen- 
sional variables. Thence the wavespeed is c = co/k and we define the wavenumber set 
K = (fc, n, c) for convenience so that K with all positive elements refers to a downstream 
travelling wave. The velocity field u is expressed as a sum of harmonic, radially varying 
traveling waves. 



where, for convenience, we define the notation K ■ x = kx + n9 — ivt (which is not a dot 
product). As such, Uk for a particular K maps a point in r to the appropriate complex 
Fourier coefficient of the velocity field at that point. 

As yet, only the wall- normal direction remains untransformed; the problem is to find 
a suitable basis for these functions of r. We introduce a projection onto a divergence free 
basis {^mir)} in the; radial direction, Fourier modes in the homogeneous spatial directions 
and in time, so that only harmonic forcing and response is considered. With coefficients 

{Xuik}, 



Re = 



V 




(2.2) 





(2.3) 
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In order to eliminate pressure, the basis functions are required to have the special prop- 
erties 

{L,^b)r = Sab, (2.4) 

V • ^{r)mK = 0. (2.5) 

The mean equation falls out naturally as the steady Fourier component of the Navier- 
Stokcs equations (the component for which for which K ■ x = always and everywhere) . 
This component of the velocity field is essentially the velocity field averaged over all 
homogeneous directions. Knowing this component implies correct Reynolds stresses in 
the rest of the model via its appearance in the resolvent. Accordingly, an assumed mean 
is used to close the equations. In principle, another Fourier component could have been 
chosen instead. It is convenient to denote this mean velocity field by uq. The field with 
the mean subtracted is denoted u = u — Uq . 

The term f = — u . Vu described of the nonlinear, triadic interaction represented 
as an internal forcing to an otherwise linear system. It is similarly decomposed into 
constituent {k and its {k,n,Lu) ~ (0,0,0) component similarly denoted fg. This finally 
yields equations for the fiuctuations that are linear in Ufc, and a base flow equation, 

-iujUK = jCKUK + iK, V(fc,n,w) ^ (0,0,0), (2.6) 

= fo - uo • Vuo + ;^V2uo. (2.7) 

The response of the flow at a particular wavenumber combination subjected to the har- 
monic forcing fff is given by rearrangement of equation 2.6, 

UK = i-ioj - CKr^fK (2.8) 

' i{kuo -u))- Re-^D -2inr~^Re-^ 1 

2inr-'^Re-^ z(fcuo - w) - Re'^D 

-drUo i{kuo-co) - Re-'^{D + r-'^) 

(2.9) 

with D = + r~^dr — r^'^{n'^ + 1) — fc^, and the states being the radial, azimuthal 
and axial velocities expressed in a divergence-free basis. The operator {—ico — Ck)~^ is 
called the resolvent and is the focus of our analysis. For large enough Reynolds number, 
we expect high gain in regions of high shear where the drUQ is large and at the critical 
layer where c is equal to the local mean velocity. This localisation is the physical basis 
for there being a highly selective response to forcing and we posit that it will be somehow 
present in any simplifying theory of coherent structure. 

2.2. Most amplified modes 

McKeon & Sharma (2010) use the Schmidt decomposition of the resolvent (the singular 
value decomposition in the discrete case) at particular K corresponding to those observed 
in real wall turbulence to rank the response to forcing in an energetic sense, i.e. 

oo 

{-iuji - Cjc)-' =Y,i^iK{r)aiK(^lK{r) (2.10) 
1=1 

with the orthogonality condition 

{(piK{r), 4>mK{r))y = 5im, {'>PiK{r), 'ipmKir))y = 5im (2.11) 
and ordering aiK > <^l+i,K > 0. 
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The inner product on r (equivalently y) only is required, 

{aK,bK)r= [ a*Kir)bK{r)rdr. (2.12) 
Jo 

This defines the forcing (0) and response (tp) modes associated with the response of the 
flow. Note that these are normaUsed with respect to the energy integrated over the radius 
of the pipe. The velocity of a mode at any given wall-normal location will therefore be 
dependent on the allocation between velocity components and the wall-normal location. 
The (piK and ipiK form the right and left Schmidt bases for the forcing and velocity 
fields at a given K and the real (Tik are the associated singular values. Assuming Ck is 
always stable, this decomposition exists for real w. The decomposition is unique up to a 
pre-multiplying unitary complex factor on both bases corresponding to a phase shift and 
up to the ordering in I of crj/cs, hence in what follows the relative phases are fixed with 
respect to the first coefficient. 

This basis pair can then be used to decompose some arbitrary forcing and the resulting 
velocity at any particular Fourier component 

oo 

fifW^^^Z-iicWx/if (2.13) 
(=1 

oo 

UK{r) = ^(TiK'4}iK{r)xiK- (2.14) 

1=1 

The energy of the same Fourier component of the resulting disturbance velocity is 

oo 

Ek = {nK, vlk) = J2 ""^kxIk- (2-15) 

1=1 

The rank-rn approximation of the resolvent at a given K by Yl'iLi ''Pik<7ik<I>*k optimal 
in the sense of the energy defined in (2.15) for equal XiK- The forcing shape that gives 
the largest energy gain at a particular K is given with xiK = for Z 7^ 1. This approach 
describes the dependence of maximum energy amplification on the radial form of the 
forcing in the wavenumber and frequency domain. The norm of the resolvent at K 
induced by the inner product over the radius is the leading singular value, (Jik- This 
means that the normalised harmonic forcing that gives the largest disturbance energy 
in the (•, •)^ sense is fK = 4>ik, with a gain of aiK- The next largest arises from {k = 
<t>2K and so on, at a particular wavenumber pair and frequency. The corresponding fiow 
response modes are given by the related u^- = ipiK, ^2K, etc. 

For the purposes of this work, we consider only real uj, k and n. This corresponds to 
statistical homogcmcity in those directions. Consideration of transient or spatially devel- 
oping fiows would require complex frequency or wavenumber respectively. We remark 
that consideration of only real ui gives a great simplification of the problem. Indeed a 
further simplification is possible: as in the streamwise direction, where projection into a 
streamwise periodic unit limits the lower range of k, the idea of a periodic decomposition 
in time implies a fundamental frequency setting the lower limit in oj. The relationship 
between the implied "turbulence period" and the unit length is given by the timescale 
associated with the largest structures captured in the box. Similarly, we need not con- 
sider u) at frequencies above those where viscosity dominates. The periodicity also results 
in discrete lo. The equivalent of this in the wall-normal direction is a natural truncation 
where the amplitude of a response mode, namely the product of the singular value and 
the magnitude of the forcing, is below a given threshold e, i.e. cTj {—ipa " VV'b, (f>i) < £• 
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As in McKeon & Sharma (2010), the computational analysis of £k was performed 
using a modified version of the spectral code of Meseguer & Trefethen (2003). The tur- 
bulent mean profiles for the higher Reynolds number were determined directly from 
experimental data from the Princeton/ONR Superpipe and reported by McKeon et al. 
(2004). Unless stated otherwise, all results shown here will be for Re = 75 x 10^ (i?+ = 
Dut/{2v) = 1800), where Ur = \/tw/p is the friction velocity, Tw is the mean wall shear 
stress, p is the density and v is the kinematic viscosity. The lower Reynolds number 
considered in Section 4.4 uses the mean profile from den Toonder & Nieuwstadt (1997) 
at Re = 10^ {R+ = 315). 

2.3. Summary of analysis 

Here, the steps in the analysis are summarised: 

(a) Fourier decomposition in symmetric directions (axial, azimuthal, time), 
(6) decompositon of nonlinear forcing and linear resolvent operator, 

(c) assumption of the zero-frequency, zero-wavenumber Fourier mode to close the 
model, which turns out to be the turbulent mean profile, 

(d) decomposition of resolvent operator by singular value decomposition to select max- 
imal gain directions consistent with the mean, 

(e) restriction to real frequencies and wavenumbers, neglecting transients, consistent 
with statistically stationary, homogeneous turbulence, 

(/) choice of cixamplc wavcmimbcrs, frequencies and mode numbers to give modes 
representative of observed structures of interest, 

(g) selection of a number of leading singular values (in the case of this paper, one) for 
each chosen response mode corresponding to the desired model fidelity, 

(ft) choice of amplitudes and phases for constructed mode combinations. 
Of these, only the first four steps are critical to the model formulation. The remaining 
steps are chosen for reasons of simplification or exposition. At the moment the choice 
of amplitudes and phases for mode combinations is the result of informed intuition. In 
principle, this choice is governed by the nonlinear interaction between modes, and work 
is ongoing to formalise this choice. 

3. Characterising the skeleton of wall turbulence 

In order to describe the structure that has been observed in turbulence, and to aid 
comparison with the literature, we examine wavenumbers and convection velocities that 
are representative of important spectral results in the literature. Most key features of 
the chosen combinations persist away from the precise values of parameters. That is, the 
structural features to be demonstrated are quite robust to the wavenumber, frequency 
and amplitude shown. Response modes at these wavenumbers form a "skeleton" of wall 
turbulence that can describe such structure in a compact way. Close to the wall, i.e. in the 
logarithmic region of the mean velocity and below, we assume at least some semblance 
of universality, such that the results from all the canonical wall flows can be interpreted 
to inform the synthesis of dominant mode combinations. In the following, the choice 
of various combinations of streamwise and spanwise wavenumbers, relative streamwise 
convection velocities and mode amplitudes is outlined and justified. 

3.1. Selection of representative wavenumber combinations 
At sufficiently high Reynolds number, wall flows are known to exhibit at least two en- 
ergetically important combinations of lengthscales, namely those associated with the 
near- wall cycle centered with energetic contribution concentrated around y+ = 15, 
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(A+,A+) « (1000,100), and the so-called Very Large Scale Motions (VLSMs, Kim & 
Adrian (1999)). The importance of the latter has been emphasized only in the past 
decade (Smits et al. 2011). The recent comparison between canonical flows performed by 
Monty et al. (2009) describes the slight differences in the VLSM lengthscales between the 
canonical flows, but for the purposes of this manuscript, we adopt a scale of \x/5 — 27r, 
such that k = 1. The spanwisc lengthscale associated with the VLSMs has been investi- 
gated by Hutchins & Marusic (2007a), Bailey & Smits (2009) and Monty et al. (2007) and 
found to be approximately equal to the outer lengthscale, n « 6. The emergence of these 
two energetic spanwise wavelengths are also consistent with the linear amplification anal- 
yses of, e.g. del Alamo & Jimenez (2006) and Cossu et al. (2009) in channel and boundary 
layer flow respectively, which reveal that the globally most amplified perturbations are 
associated with two modes with inner and outer scaling. 

Thus we expect there to be two combinations of streamwise and spanwise lengthscales 
which can be considered to constitute the spectral "skeleton" of turbulence in the near- 
wall and logarithmic region: namely (A+, A+) ~ (1000, 100) and {k, n) = (1, 6). We focus 
here on the latter combination, centered in or around the logarithmic region, but note 
that there are likely similar mechanisms to those discussed below associated with the 
near-wall motions. We will discuss the appropriate convection velocities in the following 
section. 

The recent work of Dennis & Nickels (2011a, 6) used conditional averaging of time- 
resolved, cross-stream PIV to reconstruct hairpin packets with spanwise dimension z/5 
0.5. That study used a much greater wall-normal extent than those identified in Adrian 
et al. (2000a), which was at lower Reynolds number and in pipe flow, and a streamwise 
extent of at least Q5. These packets were observed to span long low-momentum regions, 
assumed to be associated with the low-speed part of the VLSMs. A typical streamwise 
vortex spacing of approximately 5 suggests consideration of a mode combination with 
(fc, n) — (6. ±G) (mode Ka) with the VLSM. This combination is termed here the "ideal 
packet", which we represent with packet Kb- 

We will also consider a shorter, narrower mode with {k, n) — (7, ±12) such that there 
are consistent triadic interactions with {k,n) = (6, ±6) and {k,n) = (1,±6). Packet Kc 
represents what we will ultimately call the "modulating packet" , in which the modes are 
aligned in the spanwise direction. The "decorrelated modulating packet" , Kd , contains 
modes that are still triadically consistent but have a less simple azimuthal periodicity. 
Of course, in real flows, these will be subject to significant spread or jitter around these 
representative values, perhaps best interpreted in the context of a continuous spectrum 
of energetic scales. 

Wc! note also that studies that rank scales based on their contribution to the energy 
integrated in the wall-normal direction, e.g the proper orthogonal decomposition (POD) 
of Hellstrom et al. (2011), identify scales that are even longer in the streamwise direction 
than the VLSMs. Thus we also consider the extremely large-scale mode corresponding to 
Hellstrom et al. (2011), which, as will be shown below, corresponds to (fc,n) = (0.3, ±3). 
We term these "globally energetic" modes and represent them with packet Ke. 

The modes under investigation in this paper, which should be considered as represen- 
tative of the skeleton of wall turbulence, are summarized in table 1. 

3.2. Choice of convective velocity 

While a precise definition of "coherent structure" has proven elusive over the years, 
at minimum it must have the property of coherence over significant space and time. 
With this in mind, the streamwise convective velocity c for each mode was selected with 
reference to two constraints: 
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k n c A 



Ka: 


single response mode 




6 


±6 


2/3 


-l.OOi 


Kb: 


ideal packet 


Kt 
K2 


6 
1 


±6 
±6 


2/3 
2/3 


-l.OOi 
4.50i 


Kc: 


modulating packet 


ifi 
K2 
Ks 


6 
1 
7 


±6 
±6 
±12 


2/3 
2/3 
2/3 


-l.OOi 
4.50i 
0.83i 


Kd 


decorrelated modulating packet 


Kl 
Kl 


6 
1 
7 


±8 
±6 
±14 


2/3 
2/3 
2/3 


-l.OOi 
-4.50i 
0.83i 


Ke: 


globally energetic modes 


K4 

Ke 
K2 


0.3 
1.5 
2.1 
1 


±3 
±4 
±5 
±6 


2/3 
2/3 
2/3 
2/3 


1.00 
2.00 
2.00 
1.00 



Table 1. Representative mode combinations for each wavenumber packet. A is the complex 
amplitude of the response mode {Ai = aiXi as per equation 2.14). The phase angle of complex 
amplitude A determines the relative phase of that mode. 



u V w 




Figure 2. Amplitude and phase variation over wall-normal position of i!iK{r), for single re- 
sponse modes Ki (top row), K2 (middle row) and K3 (bottom row) each at unit magnitude. 
The amplitude is shown by the solid linos (lower abcissa) and the and the phase, in tt radians, 
by the dotted lines (upper abcissa). Components of streamwise, radial and azimuthal velocity 
are shown in the left, middle and right columns respectively. 
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(a) the relative velocity of one mode relative to the other must be consistent with 
observations of the evolution and dispersion of real hairpin packets; 

(b) oj must be sufficiently low for each mode such that the footprint of the mode reaches 
down to the wall, since the majority of observations of the structure under consideration 
here are have been assumed to be attached to the wall in the sense of Townsend. 

The latter constraint is simply met by examination of the velocity and swirl fields asso- 
ciated with the selected modes once the frequency has been selected. The hairpin packet 
modes recorded in table 1 all meet these criteria. 

With regard to the first constraint, we require limited dispersion of the assembled 
wavepacket, akin to enforcing similar streamwise modal convection velocities. In terms 
of hairpin packets, Adrian et al. (2000 &) reported less than 7% dispersion of packets in 
their PIV study, equivalent to noting that individual hairpin structures convected with 
similar streamwise convection velocities. This is mandated in some sense for the packet 
structure to be observable in a turbulent flow: if the difference is too large between 
convective velocities, then any packet structure observed will be highly dispersive, with 
such a short lifetime that it will be unlikely to be detected in a real flow. 

In the extreme, the validity of Taylor's hypothesis of frozen turbulence over times and 
distances which are small relative to the largest scales in the flow (equivalent to a zeroth 
order functional expansion) implies zero dispersion at a given wall-normal location and 
therefore wave packet with all components traveling at the same convective velocity. 
In order to investigate only the most persistent structure (therefore most likely to be 
observed in a real flow), we opt to follow Taylor's hypothesis in our selection of convection 
velocity c and identify structure associated with idealised, non-dispersive wavepackets. 

Figure 3 shows the variation of the wall-normal location of the peak in streamwise 
velocity intensity associated with the component modes of Kb from table 1 as convective 
velocity, c, of each mode is increased. The appropriately-scaled mean velocity profile is 
also shown to highlight the transition of each mode from a configuration in which the 
energy can be considered to be "attached" to the wall, i.e. the footprint reaches down 
to the wall and the location of the peak energy varies little with increasing c, to a 
critical mode where the wall-normal locus of the peak tracks the wall-normal location 
where c corresponds to the local mean velocity. The selection of modes with the same 
convective velocity is akin to picking modes with peak streamwise intensities that fall on 
the horizontal line in figure 3. A convection velocity corresponding to the outer (VLSM) 
mode, c = 2/3, discussed in McKeon & Sharma (2010), was selected, and appears to 
correspond exactly to the critical layer for K2 and close to the critical layer for Ki. 
Inspection of figure 3 reveals the wall-normal location of the critical layer associated 
with the modes is about y"*" = 140; it will be seen that this is a fairly good predictor of 
the location of the hairpin head. 

The effective enforcement of Taylor's hypothesis in this work has equivalence to the 
projection of use of cross-stream, time-resolved PIV into the downstream spatial dimen- 
sion using Taylor's hypothesis, for example in the work of Dennis & Nickels (2011a) 
and Hellstrom et al. (2011). In the former study, averaging conditioned on swirl criteria 
was employed, effectively winnowing away activity without time-correlation and leaving 
structure that is coherent in time. The resulting structures can be viewed as reflecting 
the effect of filtering the full flow field based on the identification of structure and char- 
acterizing the dominant structural scale distribution in fc, n and y. Use of a different 
filtering event would result in a different identified field. By contrast, the current ap- 
proach works in the opposite direction, effectively filtering on k and n and determining 
the y distribution of velocity associated with the response mode, permitting subsequent 
identification of any associated coherent structure. 
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Figure 3. Variation of the radial location of peak streamwise velocity fluctuation for the re- 
sponse modes associated with the wavenumbers of Kb (table 1) with increasing convection 

velocity, c. The lines denote — • — : (fe, n) = (6, ±6) and : (fc,n) — (1,±6). The mean profile 

at this Reynolds number, = 1800, is given by the solid curve. The horizontal line identifies 
c = 2/3, identified as the location of the critical layer for the VLSM by McKeon & Sharma 
(2010), corresponding to a wall-normal location of j/"*" w 140 at this Reynolds number. 

3.3. Mode amplitudes 

Observations of the full, three-dimensional {k, n, oo) spectrum are relatively rare in the 
literature, but the measurements of the streamwise velocity by LeHew et al. (2011) give 
some guidance on the appropriate selection of relative mode amplitudes (which will be 
important when summations of modes are considered). These data suggest that the 
amplitudes of the streamwise velocities of the turbulence skeleton should be of the same 
order of magnitude, with the highest amplitude associated with the outer mode (VLSM). 
The selected amplitudes are given in table 1. 

3.4. Tools for identification of coherent vortical structure 

The challenges associated with identifying coherent vortical structure in a turbulent 
field have been well-studied, e.g. Jeong & Hussain (1995), Chakraborty et al. (2005). 
The swirl, A, is given by the imaginary part of the complex conjugate eigenvalue pair 
associated with the velocity gradient tensor. In what follows, we will use the swirl to 
identify regions of rotation as opposed to the combined influence of shear and rotation 
identified by vorticity; note that the use of any of the other common identifiers of coherent 
structure, e.g. Chakraborty et al. (2005), give similar results in this simplified flow model. 

4. Structure from response modes 

We begin our discussion of the coherent structure associated with the response modes 
by investigating single modes, before demonstrating that structures such as hairpin pack- 
ets arise naturally from the linear superposition of such modes. 

4.1. Structure associated with individual wall modes 
The streamwise and wall-normal velocity components, u{y) and v{y) and the Reynolds 
stress, uv{y), of the response mode associated with Ka are shown in figure 4. Note, the 
distinct features of the various velocity components, namely the inclination of isocontours 
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(a) (b) (c) 




Figure 4. Heatmaps of the velocity distributions associated with Ka at R = 1800, showing 
(a) u, (b) V and (c) —uv. The azimuthal component (at this azimuthal location) is negligibly 
small. 



of the wall parallel components in the downstream direction and the corresponding lack 
of phase variation in the wall normal direction in the wall normal velocity. The aspect 
ratio of the response mode is approximately proportional to the ratio of the wall-normal 
distance to the peak intensity to half the wavelength, and this is also reflected in the 
inclination angle of the wall-parallel components to the wall. Thus different combinations 
of wavenumbers will necessarily lead to different inclination angles of coherent regions in 
the velocity field. We return to this point later in this section. 

Figure 5(a) shows the distribution of isocontours of swirl associated with the Ka 
response mode, coloured by the sense of the azimuthal vorticity. The periodic distribution 
of the velocity components give rise to a periodic array of coherent, downstream-inclined 
regions of swirl with alternating sense of azimuthal rotation in the head. There is obvious 
similarity between the geometry of these structures and the hairpin vortices observed in 
real flows via conditional averaging. Both prograde and retrograde hairpins (vortices with 
azimuthal rotation in and counter to the sense of the mean shear) respectively, have been 
identified in the literature, e.g. Wu & Christensen (2006), and similar phenomena have 
been reported by e.g. Falco (1977), Carher & Stanislas (2005). The equal distribution of 
both types here is a natural consequence of not including the turbulent mean velocity 
in the calculation, since the vorticity field associated with each mode must necessarily 
integrate to zero over a wavelength. 

A more recognizable distribution of vortices is obtained by superposing the mean ve- 
locity profile and the velocity response mode, as shown in figure 5(b). For the relative 
amplitudes of the response mode and mean profile and the swirl threshold selected for 
this figure, there is no observed signature of the retrograde vortices. While the formu- 
lation permits linear superposition of the velocity modes, the nonlinear nature of the 
swirl operator is refiected in the suppression of retrograde vortices and strengthening 
of prograde ones. Note that the apparent strength of this effect can be altered by the 
choice of different amplitude and swirl threshold combinations. McKeon et al. (2010) and 
Dennis & Nickels (2011a) have pointed out that the operations of temporal averaging 
and calculating swirl are also not commutative; swirl(ave(u(i))) ^ ave(swirl(u(t))). 

Consider the azimuthal component of two-dimensional swirl. In Cartesian coordinates 
(to simplify the discussion) this is given by 

^ 1 _^ I f du dv\^ f dudv dudv\ 

This can be used to identify hairpin vortex heads in the streamwise wall-normal plane. 
It also clearly identifies the nonlinear dependence of the diagnostic on the velocity and 
local velocity gradients. Here we examine the importance of this effect in the presence 
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Figure 5. (a-b) Isosurfaces of constant swirling strength at 50% of the absolute maximum over 
the volume for Ka (three wavelengths are shown in the streamwise and two in the spanwise 
directions), coloured with the azimuthal vorticity. Red and blue denote rotation in and counter 
to the sense of the classical hairpin vortex respectively, (a) Equal numbers of prograde and ret- 
rograde vortices when just the velocity response mode is considered, (b) With the mean velocity 
profile added, i.e. effectively under a Galilean transformation with constant (zero) convection 
velocity subtracted throughout the field of view, the retrograde vortices are suppressed and the 
prograde ones are strengthened. The ratio of center line mean velocity to is 70 : 1. 
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of a complex velocity field. Consider the swirl associated with hairpin heads as given 
by equation 4.1. For a sufficiently large contribution from the mean profile to the |^ 
term appearing in the second bracket on the right-hand side, the term under the square 
root becomes less negative, ultimately yielding only a real component and therefore no 
retrograde swirl when > (since ^ > everywhere). The opposite effect is true 
when 1^ < 0, and the apparent strength of the prograde structures is increased. Since 
1^ < through a prograde hairpin head, the swirl field under the Galilean transformation 
(effectively a subtraction of a constant, zero velocity from an equivalent real field with 
mean shear) in figure 5(b) is contaminated by the effect of the mean shear. Conversely, 
the Reynolds decomposition of figure 5(a) (the local mean velocity is absent) shows an 
even distribution of pro- and retro-grade vortices. We note that Adrian et al. (2000a) 
have discussed identifying structure using both a Reynolds decomposition (subtraction 
of the local mean velocity from an equivalent real field) and a Galilean transformation 
(subtraction of a constant, zero velocity from an equivalent real field). 

This effect is further-reaching, however, since the (linear) superposition of velocity 
response modes will lead to local velocity gradients that will alter the swirl function in 
the same way. We will exploit this effect in the reconstruction of assemblies of hairpin 
vortices into packets. 

In McKcon & Sharma (2010), we showed that the velocity response modes have charac- 
teristics consistent with real turbulent flows. The simple development above demonstrates 
that a single attached velocity response mode leads to structural features highly remi- 
niscent of hairpin vortices, while the superposition of the mean velocity profile (or use 
of a Galilean transformation) effectively reduces the observations of retrograde vortices. 
These effects are robust with respect to the selection of relative amplitudes and plotting 
thresholds. 

The inclination to the wall of the vortices associated with Kx is around 10°, consid- 
erably less than the average values of approximately 25 — 45° observed in real flows, 
however we note that the aspect ratio of the structure is a simple consequence of the 
K combination selected; therefore we proceed with Ki as the base response mode and 
build up the complexity of the wavenumber packets to observe the structural features 
associated with the superposition of multiple response modes. 

4.2. Formation of "ideal" hairpin vortex packets 

We now consider structure arising from the wavenumber-frequency combination identified 

as Kb in table 1. It will be shown that these modes, consisting of matched azimuthal 
wavenumbers, azimuthal phase and streamwise convection velocities generate "ideal", 
non-dispersive, azimuthally-aligned hairpin packets. 

Isocontours of swirl for the two-mode combination Kb are shown in figure 6. The su- 
perposition of modes that in isolation generate an array of hairpin vortices with period 
determined by K leads to a more complex pattern of coherent vorticity which is remi- 
niscent of a hairpin packet inclined in the downstream direction. In the absence of mean 
shear, packets of both prograde and retrograde vortices are observed, with a lengthscale 
determined by the longer value of streamwise wavenumber in the combination. As ex- 
pected, the effect of adding mean shear (figure 7) is to suppress the retrograde vortices, 
almost eliminating the retrograde packets for this ratio of mode amplitude to mean shear. 
Due to the ratio of streamwise wavenumbers (6:1), these packets appear to contain six 
vortices, with the last vortex in a packet occurring close to the streamwise location of 
the first vortex in the following packet with the same sense of swirl. 

Figures 8 and 9 show the velocity and swirl distributions in cuts on the plane of sym- 
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Figure 6. Isosurfaces of constant swirling strength (33% of maximum over the volume) for the 
mode combination Kb^ coloured by wall-normal vorticity. The lower two views show where the 
cuts for the next two figures are taken: the cut at constant Q for figure 8 and at constant for 
figure 9. 
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Figure 7. Isosurfaces of constant swirling strengtli (50% of maximum over the volume) for the 
mode combination Kb in the presence of mean shear, coloured by wall-normal vorticity. The 
ratio of the centerline mean velocity to |^i| is 1000:1, to \A2\ is 1000:4.5. All hairpin heads have 
prograde rotation. 



metry of the hairpin heads (at one quarter of the spanwise wavelength) and at = 45 
(through the trailing legs of the hairpins), respectively, in the absence of the mean shear. 
The distributions of in-plane velocity shown in the top panels of figures 8 and 9 for the 
streamwise/wall- normal and wall-parallel planes show that while there is some footprint 
of the shorter mode on the larger mode response in the streamwise component (filled 
contours), it is the dominant contributor to the wall-normal and azimuthal velocities 
(line contours) for this combination of response mode amplitudes. The maximum swirl 
associated with the isolated shorter mode far exceeds that of the longer mode because the 
velocity gradients are larger, even though the peak velocity amplitudes are much smaller. 
These composite velocity lobes lead to swirling motion in both planes, consistent with 
the single mode results of figures 5, but now clustered into "packets" of vortices with 
like-signed azimuthal vorticity (figure 8(b)). The main effect of the superposition of the 
large and small streamwise wavelength modes is the cause an apparent variation of the 
wall-normal location of the heads of hairpin vortices associated with the lengthscales of 
the shorter mode because of the large-scale velocity gradients associated with the longer 
one. The average diameter of the cores of the hairpin heads appears to be approximately 
50 viscous units, although this has some sensitivity to the swirl threshold selected. 

The lower panel of figure 9 gives an indication of the location of the prograde swirling 
motions relative to the low momentum regions of the VLSM mode; figure 10(a) shows 




Figure 8. (a) Heatmap of streamwise velocity for the mode combination Kb overlaid with 
isocontours of wall-normal velocity in the streamwise/wall-normal plane centered at the hairpin 
heads, n26/2-K = 0.5. (b) Corresponding two-dimensional vector plot with isocontours of swirl 
at 50% of the maximum over the slice, coloured by the sense of azimuthal vorticity. We identify 
a packet as, for example, the succession of prograde vortical motions (red) in rising formation 
from a;"*" ~ to ~ 10000. The packets of pro- and retro-grade hairpins form due to the local 
shear provided by the long streamwise streaks of K2- These act to organise the vorticity and 
modulate the shorter modes. 




Figure 9. (a) Heatmap of streamwise velocity for the mode combination Kb overlaid with iso- 
contours of azimuthal velocity in the streamwise/spanwise plane at = 45. (b) Corresponding 
two-dimensional vector plot with isocontours of swirl at 65% of the maximum over the slice, 
coloured by the sense of wall-normal vorticity. Note that the view is effectively looking from 
'underneath' figure 6. 



a more visual representation. The low speed momentum region appears to exert an 
organizing effect on the structures, such that they straddle it, tying the hairpin packet 
to occur at the same streamwise location as the low speed region. The results from 
this two-mode model closely resemble the conditional structure identified by Dennis & 
Nickels (2011a) and reproduced here in figure 10(b), although note that the condition in 
the latter data was set further from the wall than where we observe spanwise swirl for 
Kb- Recall that, in the absence of mean shear, alternating pro- and retro-grade vortices 
are associated with this K combination, so that the initial relative phase between the 
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(b) 

Figure 10. The relative position of swirl (isosurfaces in black) and isosurfaces of streamise ve- 
locity fluctuation (red and blue for high and low momentum, respectively) ; (a) model prediction 
of the swirl field (50% of maximum) and streamwise velocity (±50% of maximum) arising from 
the mode combination Kb, with mean shear of relative amplitude 500 and (b) reproduction 
of figure 18 from Dennis & Nickels (2011a) showing structure conditionally-averaged on occur- 
rence of a spanwise swirl at y/S = 0.42. While exact scale comparisons are not meaningful, the 
arrangement and mechanism appears to be similar. 

large and small streamwise scales is not the cause of this collocation in space. Rather, the 
large wall-normal extent of the streamwise velocity variation associated with the large 
scale response mode, and in particular the shear, preferentially biases the swirling 
motion to be prograde and occur along the location of maximum shear corresponding 
to the large scale low momentum region. This observation explains the relative locations 
of experimental observations of structure and large scale low momentum regions, e.g. 
Meinhart & Adrian (1995); Adrian et al. (20006); Ganapathisubramani et al. (2003), 
however it essentially reverses the causality proposed by Head & Bandyopadhyay (1981), 
who comment on shear layers consisting of a forest of inclined hairpins. Note also that 
figures 9 and 10(a) hint at the well-known result of spanwise scale growth with wall 
distance, which is at least approximately linear in the overlap region of the mean velocity, 
detailed by attached eddy theory, experimental data on the spanwise scale associated with 
conditionally-averaged eddies (Tomkins & Adrian 2005) and the spanwise velocity spectra 
of the streamwise velocity from the direct numerical simulation (DNS) of Schlatter et al. 
(2009). 

Of course, we have considered here "ideal" combinations of K, in the sense that the 
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spanwise wavenumbers and phases, and the streamwise convection velocities of the modes 
are matched, leading to azimuthal symmetry of the vortical structures and packets of 
structures that do not evolve in time. These constraints are likely the exception rather 
than the rule in real flows. However, examination of these idealized structures has permit- 
ted a demonstration of the swirling motion associated with individual near- wall response 
modes from the model and shed light on their organization duo to the interaction of 
response modes. By including a wider range of Ks, it is straightforward to obtain more 
complex swirl fields, which, not surprisingly, are harder to interpret. 

4.3. Structure self-organization, skewness and amplitude modulation of small scale 

fluctuations by the large scales 

In the preceding sections, it was demonstrated that a superposition of response modes 
could lead to clustering of vortices, while addition of the shear associated with the mean 
velocity profile weakened or suppressed retrograde vortices, leading to recognizable pack- 
ets of prograde vortices. Here we demonstrate that the latter phenomenon can arise for 
subsets of modes with appropriate phase relationships in the absence of mean shear, 
effectively giving a mechanism by which structure can self-organise, and that in such a 
scenario the mode combination also reproduces other statistical results associated with 
wall turbulence, including the basic profile of skewness of the streamwise velocity fluctu- 
ations and the apparent amplitude modulation of the small scale turbulent activity. 

Consider the "modulating packet" of case Kq in table 1. This consists, essentially, 
of the Kb modes with the addition of a further small scale wavenumber combination 
Kz = (7, ±12,2/3), chosen such that beat frequency between the smaller scale modes 
occurs at the same K as the appropriate VLSM mode, = (1,±6, 2/3). Equivalently, 
this is a consistent set of modes in the sense of nonlinear triadic interaction. All modes 
travel at the same convective velocity and we have chosen a priori a phase lag of 7r/2 for 
the VLSM relative to the envelope of the small scale activity at the critical layer, based 
on observations of real flows (Chung & McKeon 2010; Hutchins et al. 2011). 

The spatial interaction of the two smaller modes can be identified by considering the 
general case of two mode pairs with wavenumber vectors Ki and K^, vlk^^^^^^ and 
viKi^^^^'^ , differing by a wavenumber 8, that is, Ki = K — 5/2 and = K -\- S/2. 
Let the streamwise velocity component, u of the two waves have radially-dependent real 
amplitudes Ai{y) and ^3(2/) and phases (pi{y) and (^3(2/), i.e. 

u{y) = {y) + {y) 

with 

UKriy) = My) (e'^^^-^+v^(y)^ +e-'^^'-^+'^'^v)]^ 
ukM = My) (e^[-^3-a=+y3(!/)] +e-i[i<:3-a=+v3(?/)]^ 
It is straightforward to verify that 

. X A3{y)-A,iy) ^ 
"^^^ = As{y) + 

Al{y) (e''^^-''+V^(y)+^V>^ ^-z[K-x+ipi{v)+Av]^ ^^i[^-x+A,p] _^ Q-i[^-x+A,p]^ ^4 ^) 

with A(f = — V'i)/2. The second term on the right-hand side identifies a traditional 
amplitude modulation associated with the beating of Ki and namely a signal with 
streamwise wavenumber k = 0.5 and phase Atp. For the simplest case with Ai{y) = 
^3(2/), this is the only component of u. The traditional measure of amplitude modulation 
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Figure 11. Wall-normal variation of (a) amplitude of the streamwise velocity component and (b) 

phase in multiples of 2n for the three modes making up Kc- The lines denote : Ki (A; = 6), 

: K2 {k = 1), • • •: K3 (ft = 7). —x—: tps — ipi. Solid and dashed horizontal lines show the 

location of the critical layer and the region where <pi w <P3, respectively. 



takes large values twice per period of the modulating signal; in terms of the large scale 

amplitude modulation terminology currently favored in the turbulence literature, this is 
equivalent to a signal with /e = A;2 = 1 and phase Aip that takes large values once per 
period. 

The wall-normal variation of the amplitudes and phases for the three components of 
Kc are shown in figure 2. For most of the radius, AKi{y) ~ AKsiy) and, for a wide 
range of wall-normal distance around the critical layer, the two shorter modes have a 
constant phase difference. Note there is a difficulty in resolving the phase far from the 
wall where the amplitude is negligible. In addition, the VLSM mode K2 has a difference 
of phase of 37r/2 at the critical layer, giving the aforementioned phase lag of 7r/2 for the 
VLSM mode with respect to the beating envelope at the critical layer. 

The structural organization associated with Kc is shown in figures 12 and 13, respec- 
tively. The location of the planes taken through the three-dimensional field is shown in 
figure 6. The upper panels show that the beating identified for the streamwise velocity in 
equation 4.2 affects all velocity components of the small scale modes. More interestingly, 
the lower panels show that the same threshold used for wavenumber combination Kb in 
figures 8 and 9 now identifies a complex swirl distribution centered on a packet with pro- 
grade swirl that appears to grow in the downstream direction. This can be understood by 
the action of the velocity gradients (preferentially || and ||) associated with the VLSM 
mode: the large amplitude beating of the two (locally) shorter modes is spatially collo- 
cated with large positive || from the VLSM, which has the same action on the swirl as 
the mean shear, as discussed in section 4 above. Regions of retrograde swirl are observed 
upstream of and closer to the wall than the leading edge of the progradc packet, and 
upstream and further from the wall than the trailing edge. We note that this distribution 
is at least qualitatively similar to the preferred orientations of prograde vortices with 
respect to retrograde ones identified by Natrajan et al. (2007) using cross-correlation 
techniques. 

In terms of the three-dimensional swirl field, shown in figure 14, it seems that the 

VLSM exerts an organisational effect that leads to discrete packets of prograde hairpins 
distributed with streamwise and spanwise period corresponding to the beating between 
the smaller modes. Note that an increase of mean shear would serve only to strengthen 
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Figure 12. (a) Heatmap of streamwise velocity for the mode combination Kc overlaid with iso- 
contours of wall- normal velocity in the streamwise/wall-normal plane centered at the dominant 
hairpin heads, n^Q 12-k — 0.5. (b) Corresponding two-dimensional vector plot with isocontours 
of swirl at 50% of the maximum over the slice, coloured by the sense of azimuthal vorticity. The 
packet of prograde hairpins is seen as the rising array of prograde vortices (red) from a::^ ~ 
to x'^ ~ 10000. Note the weakened retrograde vortex motion at the leading and trailing ends of 
the packet. 




Figure 13. (a) Heatmap of streamwise velocity for the mode combination Kc overlaid with 
isocontours of azimuthal velocity in the streamwise/spanwise plane at = 72. (b) Correspond- 
ing two-dimensional vector plot with isocontours of swirl at 65% of the maximum over the slice, 
coloured by the sense of wall-normal vorticity. 

the prograde vortices and suppress the retrograde ones rather than change the organized 
distribution. 

A phase difference of 7r/2 between the VLSM and beating envelope at the critical layer 
serves to maximize the differentiation between prograde and retrograde vortices, however 
this choice was also informed by results in the recent literature concerning the skewness of 
the streamwise velocity fluctuation (Mathis et al. 2009, 2011) and amplitude modulation 
of the small scale fluctuations by the larger scales (Bandyopadhyay & Hussain 1984; 
Mathis et al. 2009). We now investigate the variation of these two diagnostic functions 
for Kc- 

For this combination of three response modes, an expansion of the cube of the super- 
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Figure 14. Isosurfaces of constant swirling strengtli (33% of maximum over volume) for the 
modulating packet Kc, coloured by wall- normal vorticity. The local shear is provided by the K2 
mode, which acts to organise the hairpins around regions of low-speed momentum. The lower 
two views show where the cuts for the previous two figures are taken: the cut at constant 8 for 
figure 12 and at constant j/"*" for figure 13. 



posed velocity fields shows that the skewness is given by 



Suiy) 



2A1A2A3 cos [(P2{y) + vi{y) " ^■i{y)] 
[Al + Al + Alf/^ 



(4.3) 
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Figure 15. (a) Variation of the correlation coefficient, R* -, and skewness proxy, S ■ ■ ■, with 
wall-normal distance for Kc and (b) figure 1 of Mathis et al. (2011). 



where angle brackets (•) denote averaging in the sense of a;. A proxy for the skewness 
(essentially an unsealed version) is therefore 



S = cos [Lp2{y) + ipi{y) - 933(2/)] 



(4.4) 



Note that the contributions to the mean skewness from the full velocity field will always 
be governed by triadic interactions in which the three wavenumbers sum to zero, such 
that our modulating packet can be considered as one representative of a wide range of 
triads in the real flow. 

Following Mathis et al. (2011), the Pearson correlation coefficient is defined in terms of 
the large scale signal ul and the envelope of the small-scale activity, El{us) (where the 
large and small scales can be defined using a sharp Fourier filter or Hilbert transform) 



R* = 



{ulEl{us)) 
{{uD {EL{usY)f"' 



(4.5) 



The zero of this correlation coefficient in experiment can be interpreted either as zero 
net amplitude modulation (Mathis et al. 2009) or a 7r/2 difference in phase between the 
signals (Bandyopadhyay & Hussain 1984; Chung & McKeon 2010) at this wall-normal 
distance. Recent work by Jacobi & McKeon (2012) has determined that the sense of this 
phase is to give the small scale envelope the lead in the spatial sense, which is consistent 
with the conditional averaging results of Hutchins et al. (2011). The zero crossing has 
been observed to occur at the wall normal location corresponding to both the VLSM 
peak in streamwise fluctuation, explored in more detail in McKeon & Sharma (2010), 
and the region of zero skewness of the streamwise velocity fluctuation (Mathis et al. 2009, 
2011). For packet Kc, the expression for the correlation coefficient reduces simply to 



{Al{y)) {AlUy))Y'\o^[^2{y) - A^(y)] 



R 



{{Al[y)) {AIM) 



A/2 



COS [^2 (y) + (y) - <y53 (y)] • (4.6) 



Thus the expressions for R and S are identical. 

The wall normal variations of the skewness proxy and correlation coefficient for Kc 
are shown in figure 15, together with the recent results of Mathis et al. (2011) at a similar 
Reynolds number (Rer = 2800) in a zero pressure gradient turbulent boundary layer. 
Our model predicts the trends for correlation and skewness over a wide range of wall 
normal distance, suggesting that packet Kc does indeed capture the essential structure 
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associated with real wall turbulence. Near the wall, R and S'„ display zero crossings at 
y+ « 20,45, in agreement with the data of Mathis et al. (2011), which is interesting 
particularly considering the limited number of wavenumbers present in Kc, the limited 
radial extent over which these response modes are energetic, and the different flow (pipe 
vs. boundary layer) and Reynolds number {Rcr = 1800 vs. Rcr = 2800). Since the modes 
in Kc are not energetic far beyond j/+ ~ 400, any agreement necessarily ceases towards 
the core of the pipe. Between these regions R and monotonically decrease, also in 
agreement with experiment. Figure 11(b) shows that tps — fi is constant in the range 
50 ^ 2/^ ^ 300, as denoted by the dashed lines in figure 11. For Kq, this constant is 
fixed at 37r/2. Thus the phases (fi and their rates of change with y for the two shorter 
component modes of Kc (i = 1, 3) are essentially the same over a wide range away from 
the wall, such that the variations of S{y) and R(y) are both given simply by cos 1^92, i-e. 
they both depend only on the absolute value of phase of the longest scale with k2 = 1, 
ip2, for > 60. Note that both R and S can be considered to be effectively unsealed 
versions of what would be observed in the real flow containing many more than three 
modes. However the qualitative agreement with experiment is very encouraging, and 
expands on the hereto puzzling (Mathis et al. 2009; Schlatter & Orlii 2010; Mathis et al. 
2011) connection between skewness and amplitude modulation away from the near-wall 
region. 

In simimary, the agreement between important features of Su and R and experimen- 
tal data arise from the relative phases of the modes composing Kc- Recall that the 
component wavenumbers and frequencies comprising Kc were selected from structural 
considerations, however they appear to capture the essence of two well-studied statistical 
measures surprisingly well. While the form of the skewness will always be dictated by 
wavenumber triads such as Kc, a fact that becomes quite apparent when studying modes 
decomposed into K space, this agreement suggests that the apparent amplitude modu- 
lation effect can be interpreted as a consequence of the phase relationship between the 
VLSM mode and a pair of wavenumber- frequency combinations that create a part of a 
(consistent) triadic interaction. That R{yc) and S{yc) are both zero at the VLSM critical 
layer in the real flow suggests a similar phase relationship (or at least alternatives that 
sum to zero) for the other modes present in the real flow, with any significant differences 
confined to the near-wall and core regions. Away from the critical layer, the velocity field 
is determined by the relative variations of wall-normal phase associated with the different 
modes. 

4.4. Spanwise decorrelation: Clustering of coherent structures and the largest scale 

motions 

The results above demonstrated that an appropriate phase relationship between scales 
can lead to organization of coherent structure into ideal hairpin packets associated with 
the superposition of velocity response modes with matched azimuthal wavenumbers (in 
the sense that the larger n values are integer multiples of the smallest value and there 
is no azimuthal staggering of the modes). In this subsec;tion we extend the; analysis to 
investigate structural organization in the presence of spanwise decorrelation, by removing 
this matching condition, both for the hairpin vortices and the very long wavelength 
motions. 

The mode combination Kd was chosen to explore the variation of n on the modulating 
packet by considering a triadically-consistent mode combination in which the azimuthal 
phase variation is not matched as in Kc- This decorrelation of the resultant velocity 
field impacts the appearance of coherent vortical structure in two main ways, as shown 
most clearly in the two-dimensional plots of figures 16 and 17, and the isocontours of 
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Figure 16. (a) Heatmap of streamwise velocity associated with Kd overlaid with isocontours of 
wall-normal velocity in the streamwise/wall- normal plane at the azimuthal location of maximum 
azimuthal swirl, i.e. at the hairpin heads, (b) Corresponding two-dimensional vector plot with 
isocontours of swirl at 50% of the maximum over the slice, coloured by the sense of azimuthal 
vorticity. Note the pair of retrograde vortices at a;"*" ~ 4000 and x'^ ~ 6000, just before the main 
packet begins at ~ 7000. 
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Figure 17. (a) Heatmap of streamwise velocity associated with Kd overlaid with isocontours 
of azimuthal velocity in the streamwise/spanwise plane at y"*" — 72. (b) Corresponding two-di- 
mensional vector plot with isocontours of swirl at 65% of the maximum over the slice, coloured 
by wall-normal vorticity. 



swirl shown in figure 18. Firstly, the local wall-normal gradients of streamwise velocity 
associated with K2 augment the swirl associated with the heads of prograde hairpin vor- 
tices and suppress retrograde ones, in much the same manner as in the ideal, Kc case. 
Thus the simple thresholding of swirl for the velocity response modes continues to iden- 
tify vortical structure as predominantly spanning low-momentum regions. In addition, 
the decorrelation compromises the coherence of hairpin packets; the lack of azimuthal 
matching of wavenumbers and phases leads to the asymmetric footprints of wall-normal 
vorticity shown in figure 17, indicative of the appearance of more complex structures such 
as the canes and arches, etc., identified from analysis of DNS by Robinson (1991), but 
still clustered around the large-scale low momentum regions. Thus, the spatial variation 
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Figure 18. Isosurfaces of constant swirling strength (33% of maximum over the volume) for the 
decorrelated modulating packet Kd, coloured by wall-normal vorticity. Note the pair of retro- 
grade vortices at ~ 4000 and ~ 6000, just before the main packet begins at x'^ ~ 7000. 
The lower views show where the cuts for the previous two figures are taken: the cut at constant 
9 for figure 16 and at constant for figure 17. 
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Figure 19. (a) Filtered POD data (first 4 modes, Hellstrom et al. (2011)) from a iie ~ 10"* 
turbulent pipe showing contours of spanwise velocity fluctuations, clearly picking out the VLSM 
mode as dominant in the flow at j/ = 0.2R; (b) contour plot of the spanwise velocity fluctuations 
for the mode combination Ke predicted at Re — 10* and y = 0.2R. Note that the difference in 
n of the modes is responsible for the perceived joining-up of large-scale modes and the spanwise 
decorrelation effect. 

due to the mix of wavenumbers present in a real flow, a subset of which are investigated 
here, efficiently explains a number of structural phenomena. 

The very large scale motions (VLSM), with representative wavenumbers -K^vlsm = 
(1, ±10, 2/3), were examined at length in McKeon & Sharma (2010). That paper made 
predictions of the structural form of the motions and of the variation with Reynolds 
number of the wall-normal location of the peak streamwise velocity. The VLSMs were 
considered to be of special interest because they are (according to the definition in that 
work) the largest structures that are attached to the wall. As such, they scale with both 
inner and outer characteristics. 

The temporal hot-wire measurements of Monty et al. (2007) projected into the spatial 
domain identified even longer regions of apparent coherence in the wall parallel plane, 
with meandering in the spanwise sense. Subsequently, Hellstrom et al. (2011) performed 
a snapshot proper orthogonal decomposition (POD) study on temporal velocity field 
information obtained in a pipe flow. The flow was at Reynolds number Re — 12500 
and the data were obtained using cross-stream PIV, projected into the spatial domain 
using Taylor's hypothesis. The prevailing structure regenerated from summation of the 
first four POD modes at y/R — 0.2 (reproduced here in figure 19(a)) led to favourable 
comparisons with the predictions of the VLSM mode shapes in McKeon & Sharma (2010), 
with the caveat that the experimentally-measured observations seemed to be dominated 
by longer scales than the fc = 1 mode used by McKeon & Sharma (2010), even once the 
meandering associated with the superposition of modes was accounted for. We would add 
that at this wall-normal location, the VLSM and the other modes in Ke will be converting 
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at a different speed (c = 2/3) than any locally critical modes which presumably would 
be observed. Thus, care should be taken when interpreting measurements using Taylor's 
hypothesis for structure spanning a wide range of scales. 

The meandering effect described by Monty et al. (2007) has been explained in various 
ways (Sharma & McKeon 2011; Dennis & Nickels 2011a; Hellstrom et al. 2011) to arise 
from a decorrelation of the energetically-dominant VLSMs due to motions at other scales. 
Figure 19(b) shows a wall-parallel cut at y/R = 0.2 of the streamwise velocity field 
arising from the superposition of the velocity modes comprising packet Ke in table 1. 
The three modes therein were selected by eye to match the key large scale features of 
Hellstrom et al.'s results shown in figure 19(a). The decorrelation associated with the 
range of azimuthal wavenumbers, n, present in Ke naturally gives rise to the apparent 
meandering effect and also creates a spatial localisation of coherence in the azimuthal 
direction because of the azimuthal analogue of the beat effect described in section 4.3. 

4.5. Shear stress variation and comparison with models 

We now consider briefly the contributions to the shear stresses arising from individual 
response modes and the correspondence between the location of the maxima and obser- 
vations of structure discussed above. Consideration of the mode shapes in figure 4(a-c) 
emphasizes that the shear stress arising from a particular mode will be strongly depen- 
dent on the relative shapes of the streamwise and wall-normal responses. The typical 
shapes for the modes considered here involve isocontours of streamwise velocity that 
are heavily inclined in the downstream direction, i.e. have decreasing phase, (p{y), with 
increasing y, but have little or no phase variation for the wall-normal velocity. 

Consider the shear stress from an individual response mode with streamwise and wall- 
normal velocity fields with magnitudes ujciy) and vjiiy) and phases ^Pu{u) and (fy{y), 

u{y) = UK{y) (e'^^-^+'PM] + ^-ilK-a^+vM]^ (4_7) 
v{y) = VK{y) (e'[-^""+'^''(^)l + e-'['^-»=+V"(?')]) . (4.8) 
The contribution to the shear stress can be decomposed into mean and fluctuating parts, 

uv{y) = UK{y)VK{y) (ei[K-'^+'Pu{v)] Q-i[K-x+v>u(y)]^ (^^[K.x+v>.(y)] ^-i[K.x+v.{y)]^ 

= UK{y)vK{y) (e^^'p^^y'i-'p-^y^^ + c-'^vM-fviv)]^ 

+UK{y)vK{y) (^e'm-x+v^(y)+^^(y)] ^ ^-i[2K-x+^^(y)+^^(y)]^ 

SO 

uv{y) = 2uK{y)vK{y) [cos {(pu{y) - ^v{y)) + cos {2K ■ X + Lpu{y) + Vv{y))] , (4.9) 

or a contribution to the mean shear stress whose magnitude depends on the relative 
phases of u and v plus a fluctuating component at 2K. Recall that response modes give 
the velocity relaMve to the local m,ean. The contribiition to the instantaneous shear stress 
arising from a particular mode will be localized around the wall-normal location(s) at 
which u and v have a relative phase that is close to an integer multiple of tt and non- 
negligible amplitudes. For the modes of this study, these two conditions tend to be met in 
the interior of the flow, in a relatively narrow region away from the wall. This causes the 
contribution to the mean shear stress to have the type of profile shown in figure 20, which 
shows the wall normal distributions of shear stress associated with the modes comprising 
packet Kc- Note the large magnitude of shear stress associated with the VLSM mode, 
in agreement with the observations of Guala et al. (2006) and Wu et al. (2012). 
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Figure 20. Wall-normal profiles of Reynolds stress normalised with absolute value of the 
maximum associated with the following modes identified in table 1: (a) K\ = (6,6,2/3), (b) 
K2 = (1,6,2/3) and (c) = (0.3,3,2/3). 

Further inspection reveals that for the majority of response modes, the shear stress 
peak occurs in the region of the critical layer, whore u and v arc tt out of phase such that 
the peak contribution to the mean shear stress is always negative, i.e. {—uv)^^^^. > 0, and 
at the streamwise locations corresponding to the (co-located) minima in the streamwise 
gradient of v and wall-normal gradient of u. This location was also previously identified as 
the location of the hairpin vortex heads. The gradient of the Reynolds stress is therefore 
high and positive just below the critical layer and high and negative just above. We note 
that laminar, inviscid theory predicts the Reynolds stress gradient to be infinite at the 
critical layer, and zero elsewhere (locally) around it (Drazin & Reid 2004). 

We conclude the presentation of the results by noting that the connection between the 
present discussion on Reynolds stress and the structural, statistical and spectral results 
discussed above and information in real flows is greatly facilitated by the ability to linearly 
superpose response modes to obtain trends applicable to the full field. Unfortunately 
this technique is not amenable to the quadrant and VITA analyses performed by many 
previous workers, since those approaches are conditioned on events in the physical domain 
in which all modes participate and periodicity is returned only in the averaged sense. 

5. Discussion 

The extensive foregoing analysis, together with the results of McKeon & Sharma 
(2010), has demonstrated that simple linear combinations of velocity response modes 
reproduce key statistical and structural features of wall turbulence. Further, they can 
give insight into the underlying mechanisms for extraction of energy from the mean flow 
and for structural dynamics. We now briefly discuss and connect the key results, outline 
the mechanisms responsible for this success and make connections to observations and 
phenomenological models in the literature. 

5.1. Connections to the attached eddy hypothesis, edge states and experimental 

observations of structure 

The Reynolds shear stress distributions associated with linear superpositions of velocity 
response modes and the associated observations of structure are entirely consistent with 
expectations for the behavior of the Reynolds stress in fully turbulent flow. They are also 
highly reminiscent of the form predicted in the literature by the attached eddy theory 
of Townsend (1956, 1976) and Perry & Chong (1982) and subsequent authors, as well as 
the scale hierarchies of the mean momentum balance approach of Klewicki et al. (2007). 
We explore the similarities in more detail below. 

The attached eddy and mean momentum balance analyses incorporate and predict, 
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respectively, self-similax families, or hierarchies, of coherent structures that have a shear 
stress signature that is localized away from the wall. We have shown that the first velocity 
response modes for K values that arc attached to the wall (in the sense of Townsend, 
with a footprint that reaches down to the wall) constitute such structure. With respect to 
the mean momentum balance approach, this is perhaps no surprise since both analyses 
deal with the Navier-Stokes equations themselves. The velocity response modes likely 
identify the velocity and vortical structure giving rise to self-similar local distributions 
of Reynolds stress associated with the hierarchies. 

Similarly, the velocity response modes appear to provide the full velocity field as- 
sociated with the distribution of self-similar attached eddy structures of the attached 
eddy hypothesis, or the tpu function of Perry et al. (1986). Whereas the attached eddy 
hypothesis makes a phenomenological argument for the distribution of attached eddies 
that explains the statistical behavior of the velocity fluctuations near the wall, our ap- 
proach identifies structure from the most amplified response modes associated with the 
Navier-Stokes equations (incorporating appropriate boundary conditions), but lacks the 
appropriate distribution of amplitudes. This is the opposite interpretation to the vortex 
induction arguments given by Adrian (2007), the apparent disparity arising because of 
the assumption of the mean velocity profile as opposed to the attempt to determine the 
mean velocity profile in the latter case. As such, the approaches appear to be comple- 
mentary. 

The VLSM mode, K2, clearly belongs to a class of mode that, for high enough Reynolds 
number, should be considered both "attached" to the wall and "inactive" in the sense 
of Townsend (1956): the response has a velocity footprint that reaches down to the wall 
but a mean shear stress that is localized sufficiently far from the wall that it has no 
contribution close to the wall, but can contribute significantly to the overall shear stress. 
This last effect is expected to increase with Reynolds number, in agreement with the 
observations of Adrian (2007), who summarizes the large contributions to the mean shear 
stress arising from the very large scales in pipe and channel flows (from the original work 
of Guala et al. (2006) and Balakumar & Adrian (2007), respectively), and of Marusic 
et al. (2010) concerning the increasing significance of turbulence production in log region 
with increasing Reynolds number. 

With respect to the origin of the VLSMs, Kim & Adrian (1999) proposed that they 
consist of long packets of coherently aligned hairpins. The picture developed here, how- 
ever, suggests that the packet nature is likely due to a "beat" effect between modes 
with the specific beat frequency linked to the dominant energetic VLSM. Furthermore, 
although the VLSM acts to organise the hairpins, the two structures may be predicted 
separately. It is however our suspicion that the VLSM and hairpin response mode com- 
binations interact in a mutually supportive manner. The recent papers on hairpin-like 
exact solutions mentioned in the introduction, in particular that by Cherubini et al. 
(2011) in which a self-sustaining edge state looking like a hairpin packet was found, is 
suggestive that this is the case. This indicates that the roll-streak interaction mechanism 
is not uniquely responsible for driving turbulence and that an analogous hairpin-streak 
interaction mechanism may exist alongside the near-wall cycle. 

The Reynolds stress distributions associated with individual velocity response modes, 
which can be superposed, also illuminate the underlying structure behind recent results 
concerning the wavenumber decomposition of the gradient of the Reynolds stress (Guala 
et al. 2006; Wu et al. 2012). Figure 21 includes a reproduction of the net force spectra 
of figure 17 in Wu et al. (2012), consisting of the spectral decomposition of the sum of 
the gradient of Reynolds stress and a contribution from the Reynolds stress itself arising 
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the dropped dotted lines show uq/Ucl = 1/2 and uq/Ucl ~ 2/3. 



from the cylindrical coordinate system, 

, , d (uv) (y) (uv) (y) 

net force y = ^—^+^—^. 5.1 

dy R-y 

Note that in the framework of Wu et al. (2012), the full four-dimensional velocity field 
has been averaged over n and uj. In figure 21, the net force associated with streamwise 
wavenumbers representative of the VLSM modes, \x ~ 10, has a wall-normal extent 
of positive sign (marked D in figure 21) that reaches much further from the wall than 
observed for all higher wavenumbers. The change in sign for these wavenumbers occurs at 
or close to the location where Uq/Ucl = 2/3, which is the location of the critical layer for 
the VLSM, as predicted by the theory. These results indicate that such observations are a 
consequence of the VLSMs being critical modes with critical layers that exist further from 
the wall than for the shorter modes, as shown in figure 20(b), such that the associated 
region of positive gradient extends well into the overlap region of the mean velocity and 
the sign change resides at or near the critical layer. This effect dominates over the other 
nearby wavenumber combinations with different wall-normal extents simply because of 
the energetic dominance of the VLSMs. 

Some other observations regarding figure 21 are worthy of further comment. Firstly, 



34 



A. S. Sharma and B. J. McKeon 



the locations of the other obvious changes in the sign of the Reynolds stress gradient, 
corresponding to the transitions from regions B-E and A-C, are coincident with the wall- 
normal locations of mean velocities equal to [/+ « 12—15, respectively. These transitions 
correspond to the lowest convection velocity of structures observed for the near-wall 
region in the literature (a summary of the relevant literature is given in LeHew et al. 
(2011)), suggesting that this region is dominated by critical modes traveling at this range 
of minimum velocities. The B-E transition appears to also occur where Uq/Ucl = 2/3, 
identified by McKeon & Sharma (2010) as an important phase velocity associated with 
the VLSMs but this is probably an effect of the low Reynolds number associated with the 
DNS. Remarkably, certain modes close to the wall (not shown) can give rise to positive 
Reynolds stress, or positive then negative wall-normal gradient of {uv), in the same 
wavelength range identified by Wu et al. (2012) denoted F in figure 21. This is another 
topic worthy of further study. 

5.2. Structure arising from the Navier-Stokes equations and links to the classical 

laminar, inviscid theory 

We believe this work to constitute the first demonstration of the development of hairpin 
vortex packets based on analysis of the Navier-Stokes equations and mean profile alone. 
While many researchers have identified the non-normality of the linearized Navier-Stokes 
operator as a source of finite-horizon energy amplification, or transient energy growth, 
the formulation of McKeon & Sharma (2010) permits further identification of the loss of 
symmetry of the resolvent associated with the presence of the wall as the root cause of 
coherent structure. The key relevant insight is that the large response near the critical 
layer can generate localised Reynolds stresses that can in turn excite the linear critical 
layer amplification. It seems likely, therefore, that underlying all of the nonlinear solu- 
tions reported in the literature must be some combination of two linear amplification 
mechanisms: the critical layer amplification effect manifested as a near-singular response 
of the linearised system to forcing, and the interaction with the mean shear manifested 
as operator non-normality. 

Regarding non-normality, an operator X is normal when XX* = X*X and non-normal 
otherwise. Since the adjoint X* is defined only with respect to some inner product (such 
that {a,Xb) = {X*a,b)), non-normality is also defined with respect to the same inner 
product. Self-adjoint operators are therefore normal. As is well known, the pipe flow 
equations are invariant under the translation in t, x and 9. Consequently, the resolvent 
is self-adjoint under the integrals over these two coordinates. This is readily proved by 
considering the translation operator in each coordinate; it commutes with the resolvent 
and its eigenvectors are the associated Fourier modes. This results in an orthogonality 
condition and the natural use of the Fourier decomposition. The nonlinear interaction 
between Fourier modes then obeys the rules of triadic interaction. In the symmetric 
directions, therefore, the forcing and response modes are equal and are the Fourier modes. 
The presence of the wall changes this state of affairs: since the symmetry is lost in 
the direction normal to the wall, the; resolvent is not self-adjoint with respec;t to the 
inner product associated with integration over that direction. In this case, we see that 
< {4'a,4'b)r ^ 1, we lose the orthogonality property in the wall- normal direction and 
associated rules of triadic interaction, the forcing and response modes are no longer equal 
and the fluctuations may now gain energy from the interaction with the mean flow. The 
potential for momentum production due to this interaction is quantifled by this loss of 
orthogonality. 

It appears from figure 3 that the hairpin head is located near the critical layer. This is 
strongly suggestive that the hairpin heads arise from a mechanism essentially similar to 



On coherent structure in wall turbulence 



35 




Figure 22. Kelvin's "cats' eyes" streamlines associated with inviscid potential flow in the 
region of a critical layer (Drazin & Reid 2004). 

that described by the famous cats' eyes diagram of Lord Kelvin for the laminar, inviscid 

case - although the laminar, inviscid theory predicts the Reynolds stress gradient to 
be infinite at the critical layer, and zero elsewhere (Drazin & Reid 2004). The visual 
resemblance is demonstrated in figure 22 by a sketch of the corresponding local streamline 
pattern in the presence of the mean flow. In the present case, the response modes are 
three-dimensional, rather than the two-dimensional form implied by Squires' theorem. A 
hairpin head can therefore be interpreted as a pair of cat's eye vortices meeting obliquely, 
with regularization at the critical layer due to viscosity. 

In a further interesting link to the inviscid laminar theory, for k = 1 and at high Re, 
the point c = 2/3 corresponds (on the complex plane) to region where the S-family of 
eigenvalues of the laminar Ck impinges on the iuj axis. This explains its association 
with a high leading singular value. In addition, this is the region where the S-, A- and 
P-families of eigenvalues meet, and is associated with high sensitivity to perturbation 
(Reddy et al. 1993). Since at this point neither of the branches associated with either 
the centre eigenmodes or the wall eigenmodes is dominant, the response of the flow has 
characteristics of both. 

5.3. Taylor's hypothesis, convectiue uelocity, critical layer localisation and coherence 

The response modes under investigation are propagating in the streamwise and spanwise 
directions and distributed in the wall-normal direction. Their wave-like nature implies 
that the streamwise propagation velocity of each mode is given by — w/fc = c. Note 
that uq ~ c modes will dominate in any region with low viscosity. The result of this 
is that the energy of the critical-type response modes is localised in the region around 
the critical layer. This is less true for globally-energetic modes, which are less localised 
and cannot detach. Such a localisation is manifested by the narrowness of the ridge in 
the joint (k, u>) spectrum, where a narrower ridge implies lower dispersion and higher 
coherence. This coherence, when coupled with a localisation of energy in the wall-normal 
direction near the critical layer, is equivalent to an invariance when travelling in the 
reference frame moving at c. The localisation in the wall-normal direction and the low 
dispersion of corresponding structure is predicted by the radial dependence of the forcing 
and response modes. In turn, this is explained by the high response of the resolvent in the 
region of the critical layer. This effectively justifies the use of Taylor's hypothesis when 
searching for coherent structure: the application of Taylor's hypothesis is synonymous 
with coherent, critical response modes. In this sense, the model provides a predictive 
explanation of coherent striicture observed in wall turbulence. 

We have only shown "perfect" hairpin packets here, but the effects described above 
will be even more important in the case of modes with different convection velocities, 
where packets will appear to grow and decay, or simply "evolve" , as a consequence of the 
interference between modes. The apparent evolution of collections, or packets, of hairpin 
vortices therefore occurs because of relative motion between modes. A packet will grow 
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(decay) in space and time if the phase velocity associated with the smaller modes is larger 
(smaller) than that of the VLSM. Thus a spread in the distribution of convective velocities 
in the assembly of velocity response modes can capture the experimental observations 
of the evolution of vortical structure. This is fully within the scope of our approach but 
is not explored here. In real flows, deviations from perfect matching of c are of course 
possible. While a slight dispersion negates the problem of choosing an appropriate phase 
relationship between the modes, it is possible that the nonlinear feedback may lock the 
relative phase between modes. 

5.4. Self-organisation, phase and nonlinear feedback 

Turning now to coherent structure, we have shown that the velocity gradients associated 
with response mode combinations result in very complex swirl fields and that interference 
between modes results in the apparent disappearance of periodicity. 

The presence of an energetic very large scale mode such as K2 generates local wall- 
normal gradients of the streamwise velocity. This locally increases or decreases shear, 
eS'ectively increasing or decreasing the local Reynolds number. Because the swirl diag- 
nostic is nonlinear, this also appears to align hairpins (and their remnants such as canes 
and arches) to the low-momentum regions, in packets with the period of the VLSM. The 
large scale separation of the two types of structure (their separation in wavenumber- 
frcquency space) justifies similar quasi-static explanations. Further, the existence of ret- 
rograde vortices alongside the prograde ones is entirely consistent with observations in 
the literature, including the preferred orientations of prograde vortices with respect to 
retrograde ones discussed by Natrajan et al. (2007). 

The analysis also suggests the nonlinear travelling wave solutions reported in the litera- 
ture for transitional flow continue to be important in high Reynolds number turbulence. 
We have supported this idea by applying a simple method which examines the linear 
amplification mechanism that selects such structure. We find that if the forcing resulting 
from a mode combination produces, via the resolvent, a component of the original mode 
combination, then the mode must be self-sustaining above a certain amplitude. Such 
a combination would have to be triadically compatible and in the radial direction, the 
generated forcing must be at the same locality as the mode combination. 

We may go on to speculate that another role of the VLSM in a real flow is to cause 
an effective phase-locking of a wide range of the smaller scale modes. This effect would 
become increasingly dominant with increasing Reynolds number as the VLSM becomes 
more and more energetic relative to other individual modes in the flow. Indeed, it seems 
plaiisiblc that the quadratic nonlinear feedback acts to provide phase-locked loop be- 
haviour (analogous to how the nonlinear term acts to keep a soliton localised in the 
direction of travel), because the nonlinearity is multiplicative between frequencies and so 
has the right form to act as a phase detector. In this sense the nonlinear feedback can be 
considered as another mechanism by which the flow can self-organize in a relative phase 
sense. Further investigation of this phenomenon is outside the scope of this manuscript, 
but forms the basis for future work. 

5.5. A new turbulence kernel 
The preceding sections have demonstrated that the linear combination of velocity re- 
sponse modes constituting Kc captures both coherent vortical structure and the essence 
of higher order statistical results. In particular, the structure shown in figure 14 is 
highly reminiscent of the dominant organization to be expected from the hairpin packet 
paradigm that has dominated the literature in recent years, while figure 15(a) shows a 
distinct similarity between the skewness of the streamwise velocity fluctuations and the 
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apparent amplitude modulation of the small scales by the VLSMs associated with Kc 
and observations in experiment. 

The statistical picture that emerges from Kc, which was synthesized using insight from 
experimental observations at high Reynolds number, confirms that this simple kernel 
captures key aspects of the turbulence. In particular, it appears that the statistics of 
the log-law region are essentially controlled by the VLSM. The ability of the three- 
mode model to identify the similarity between the amplitude modulation correlation 
coefficient and the skewness is a consequence of the interpretation of both these measures 
(and the odd order moments) as providing information on the relative phase between 
modes, as alluded to by Bandyopadhyay & Hussain (1984) and Chung & McKcon (2010). 
This provides a missing piece of the previous work of McKeon & Sharma (2010). The 
importance of the VLSM to the even moments and the amplitude modulation has been 
well-documented by Mathis et al. (2009) for example. In the case of an energetically 
dominant VLSM, we expect triadic interactions of the sort described by Kq to dominate 
the skewness and the amplitude modulation, thus explaining the puzzle of the connection 
between them that has been discussed by Schlatter & Orlii (2010) and Mathis et al. 
(2011). These effects arise from a linear combination of modes with different phases, 
rather than a nonlinear effect. 

That the type of triadic interaction capable of capturing the skewness profile also gen- 
erates recognizable coherent vortical structure supports our hypothesis of its importance 
to wall turbulence. We suggest, then, that a subset of three modes that includes the 
VLSM and two smaller modes constitutes a new kernel of turbulence in the sense of 
Smith et al. (1991). Such kernels form studiable sub-units that contain important dy- 
namics associated with the flow, a well-known example being the autonomous near-wall 
cycle (Waleffe 1997) , and which can be summed together to regenerate the full picture of 
wall turbulence (or at least create a synthetic proxy for the real flow) . It is also likely, but 
is not proven here, that such a kernel self-sustains, so that the; forcing arising from any 
two components of triadic interactions like Kc is sufficient to excite the third component. 

6. Conclusion 

McKeon & Sharma (2010) developed a formulation of the Navier-Stokes equations 
where the resolvent acts as a directional amplifier to an imstructured forcing generated 
by the nonlinear interaction between modes. It was demonstrated that this predicts a 
velocity response mode associated with the first singular value of the resolvent at any 
given scale and that these modes capture statistical features associated with pipe flow 
turbulence. This was essentially a demonstration that resolvent is well approximated 
by a low-rank operator. This reflects the fact that the underlying dynamics are well 
approximated by low-dimensional dynamics, which can provide a vast reduction in the 
dimension of the problem. In the present work, we have built on this formulation to 
demonstrate that the response modes efficiently encode information on coherent structure 
as defined phenomenologically in the; introduction. While the response modes can be 
superposed to assemble more complex flow fields, any measure of swirling structure is 
inherently nonlinear and therefore the perceived structure associated with an assembly 
of modes is not the linear sum of the structure associated with individual modes. This, 
perhaps, lies at the heart of the difficulty of observation and characterization of coherent 
structure in real flows. 

The response modes contain coherent structures as objects mathematically derived 
from the Navier-Stokes equations under the assumption of a mean turbulent profile. As 
such, we posit that they form a desirable basis for a low-order model when compared to 
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approaches that repUcate coherent structures found in experimental data (such as POD) 
or describe them phenomenologically. The minimal-order model is derived from the lowest 
rank (rank-1) approximation of the resolvent; assuming the resolvent is compact, there 
is finite approximation error associated with this truncation. Though the inclusion of 
additional response modes may improve the approximation, this has not proved to be 
necessary for the present purpose. 

Whilst the mechanisms underlying turbulence are captured by the resolvent, to actually 
observe hairpins as structure in an experimental flow, four ingredients are necessary: (1) 
a nonlinear diagnostic such as swirl, (2) an averaging (or equivalently a frequency domain 
analysis), (3) coherence, provided by the assumption of Taylor's hypothesis or a critical 
layer response and (4) a flow boundary, to provide characteristic scales and a shear, 
which acts in the mean sense to provide net vorticity responsible for the hairpin head 
and locally in response to the VLSM to organise hairpins around areas of low momentum. 

The directional amplification of the resolvent provides a mechanism to understand the 
interaction of simple external forcing strategies, such as the introduction of a single K 
disturbance at the wall using dynamic roughness explored by Jacobi & McKeon (2011), 
as well as the good approximations to the mean velocity profile obtained using stochastic 
forcing and a two-dimensional, three-component flow model (Gayme et al. 2010; Bour- 
guignon & McKeon 2011). Together with the nonlinear feedback term, the directional 
amplification is understood to underlie the robustness of the flow features observed in 
real flows. There are obvious related implications for flow estimation and control. 

Finally, we have identified a new turbulence kernel, composed of these response modes, 
that explains some current controversy in the literature and, importantly, links statistical 
and structural observations by consideration of the relative phase between modes. 

Our study here has becm focused on i?+ = 1800, at the upper end of the range of 
Reynolds numbers reachable rising current DNS, but we stress that the only constraint 
on studying higher Reynolds numbers lies in the numerical precision required to deal 
with the large response associated with the most amplified forcing. 
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